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Abstract 

We are reporting a new non-Fermi liquid type normal phase that has a 
well defined Fermi energy, but without showing any non-regularity in the 
momentum distribution function n& in the whole momentum space, the sharp 
Fermi momentum concept being undefinable. The system contains a natural 
built in gap that is visible in the physical properties of the system at nonzero 
temperatures. The presence of a flat band in multi-band interacting Fermi 
systems with more than half filling is the key feature leading to such a ground 
state, which is not restricted to one spatial dimension and emerges in the 
proximity of an insulating phase. 
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Our understanding of the behavior of interacting fermionic many-body systems is inti- 
mately connected to the concept of Fermi-liquid introduced by Landau many decades ago 
|jlj] . In a normal state that preserves all symmetry properties of the high temperature phase 
of a fermionic system, the Fermi liquid behavior has been clearly observed in the normal 
state of He3 and simple metals |2j. It has the meaning that in spite of the inter-particle 
interactions, the low energy behavior of the system can be well described within a picture of 
weakly interacting quasi-particles ||. This picture can be also mathematically formulated 
d]. In these terms, in a normal Fermi liquid: a) There is a one-to-one correspondence be- 
tween the non-interacting one particle states and interacting single-particle states. This is 
concretely obtained by describing the interacting system using a perturbation theory that 
is convergent up to infinite order, b) The single particle Green-functions have a quasi- 
particle pole that gives rise to a discontinuity of the momentum distribution function 
at the Fermi surface whose position is specified by a sharp Fermi momentum value kp. c) 
The residual quasi-particle interactions can be described by a small number of parameters, 
called Landau parameters, which can be deduced from a microscopic theory taking into 
account non-divergent two-particle vertex functions |J. In the last decade however non- 
Fermi liquid behavior has been observed experimentally in the normal phase of a variety of 
materials, including higher than one dimensional systems of large interest. Examples are: 
high temperature superconductors 0, heavy- fermions 0, layered systems |§, quasi-one di- 
mensional conductors, doped semiconductors, systems with impurities, materials presenting 
proximity to metal-insulator transition ||, etc. These results are often discussed in terms of 
multi-band models ||, the presence of a some kind of gap in the normal phase being clearly 
established in many cases and subject of intensive studies PJlQ|. The last decade witnessed 



a huge intellectual effort |TT| for the understanding of the non-Fermi liquid behavior in the 
normal phase |12j of fermionic systems. On the theoretical side however, for pure systems, 
the existence of a non-Fermi liquid in a normal phase has been exactly proved only for the 
one dimensional case (i.e. Luttinger liquid JIB] ). So far, the possibility of extending the 
proof to two spatial dimensions has not been demonstrated rigorously. In fact, a rigorous 



theory of a non-Fermi liquid normal state in higher than one spatial dimensions is missing. 
For this reason, the theoretical understanding of different phenomena observed in the mate- 
rials listed above is relatively poor and theoretical advance in this subject is badly needed. 
Our work on the periodic Anderson model (PAM) at nonzero and finite on-site Coulomb 
repulsion (U) (a prototype of two-band systems containing strong correlation effects) was 
motivated by this state of affairs. We are presenting in this Letter for the first time, an 
exact solution for this model. The solution present in a restricted (but continuous and infi- 
nite domain of the phase diagram) represents a new type of non-Fermi liquid behavior in a 
normal phase for a system that has a built in gap. The obtained ground state energy cannot 
be expressed as a sum of contributions of the on-site Hamiltonian terms, and the ground 



state expectation value of the kinetic energy terms is nonzero and negative fL4| . The state 
emerges in the vicinity of a Mott insulating phase in a continuous domain of concentration 
above 3/4 filling and has a well defined Fermi energy, but the Fermi momentum cannot be 
defined, being without any non-regularity in momentum space. The property is due to 
the emergence of a flat band in a multi-band system with more than half filling and can 



be extended in an exact manner to two spatial dimensions ]T5|. Such features have been 



observed experimentally, e.g., in ARPES data, which, even for high-T c materials often reflect 
main bands without any sharp characteristics in nk |16| or necessitate the assumption of 



the presence of flat bands ||17|| . Band structure calculations for layered systems often show 
a Fermi level positioned exactly at the bottom of a conduction band with a relatively large 
effective mass around its minimum, below which a gap is present ||. Connections between 
the emergence of superconductivity and flat dispersions were also clearly pointed out in Ref. 
||18|| . Flat-band features are present in heavy- fermion systems as well |19j, and can be even 
produced by squashing carbon nanotubes EC 



We consider two bands denoted by b = c, /, the starting ID Hamiltonian being H = 
Ho + U. The Hubbard term isU = UJ2i MtMi an d we have Ho = T c +Tf + Ef + Hh with the 
kinetic energies T b = t b b\ a b i+lta + h.c. , the on-site / level energy Ef = Ef 
hybridization energy Hh = Vq + Vi, and 



Vq = E [voA*f* + h - c - 

Vi = E [Vi (i/i+i. + + h.c] . (1) 

We are presenting the case when H h contains imaginary coupling constants, i.e. V = i V 
and V\ = i V\, where Vo, V\ are real. For the starting point we consider 3/4 filling. In this 
Letter we shall give an exact solution of this problem, valid in a restricted domain of the 
phase diagram relevant for our study. In order to have a clear image of the solution and its 
physical meaning, let us consider the U = case first. To this end, we would like to express 
H through a, B = B l,a B i,o term , where B it<J = a c ij(T + ft c i+1>a + 7 f i)ff + 5 / i+lj(T and 
a, ft, 7, 5 are constants to be determined. Introducing the total particle number operator 
N = J2i,a,b ^Aai we fi 11 ^ then that ^0 = — B + r] N holds, if the following relations are 
satisfied 

V = \a\ 2 + \ft\ 2 , (2) 

E f = V ~ (M 2 + \6\ 2 ) , (3) 
-t c = a* ft = ft* a, (4) 

- tf = 7*5 = 5*7, (5) 
% V = a 7* + 6* ft = - ( a* 7 + 6 ft* ) , (6) 
i V 1 = a 5* = 7 ft* = - a* 5 = - 7* ft . (7) 

Eqs.([| - |7|) determine through a, ft, 7, 5 the value of Ef and 77 for which the presented 



structure of H is valid. With 77 = 2 yt> 2 + t\, v = Vq V\ / ( 2 tf) we obtain 

= 77 ( 1 - t?^r 2 ), V^ 2 = -t c ^. (8) 

Taking now into consideration that B\ a B ia + B i a B\ a = \a\ 2 + \ft\ 2 + \-y\ 2 + \6\ 2 and 
imposing the condition ( N ) = 3 L where L represents the number of lattice sites (3/4 
filling), we find f/o — E g + P , where E g = Lr\ ( 1 — 2 \m\ 2 ) and m = tf / V\ has been 
introduced. The operator P = J2i,a B\ being positive semidefinite, E g is the ground 
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state energy, and the ground state wave function is that \ip g ), for which we have P\ip g ) = 0. 
We demonstrate now that 

ivv> = n f ii F i a) ) io>, (9) 

i=l \ a=l / 

where | ) is the bare vacuum with no fermions present, D i r7 = a* ( cj CT + m f\ a J + 

l(T (a)=t,4.- In order to see this, one can easily 
verify that B^ a = — F^ Bj a independent of the indices, and B\ D i cr = 0, so 

P\ij)g) = 0, and | ijjg ) is the ground state. We must stress, that for U = the ground 
state is entirely given by f} and F[ , the F^ operator being completely arbitrary apart 
from the requirement that it introduces under the Yl { product L electrons into the system. 
The concrete expression of F{ is fixed by the nonzero U as follows. When U ^ 0, the 
Hamiltonian contains, besides -£/o> the Hubbard U as well. However, one may observe that 
U can be exactly transformed as U = U P' + U J2i,a Ma — U L , where P' = J2i P[ an d 
P[ — ( 1 — n{+ — n{\ + n{\). But, P[ is one if on the i site there are no / electrons, 
and is zero, if on site i there is at least one / electron. As a consequence, P' gives a sum 
of non-negative numbers, so it is a positive semidefinite operator. Furthermore, P 1 gives its 
minimum eigenvalue (i.e. zero) for a wave function that contains at least one / electron on 
every site of the lattice. Let us consider for this reason 

^ (3) = E a* fl . (10) 

(7 

where the a a are constants to be determined. Now | ip g ) introduces 3 electrons L times 
within the system, so the solution (two bands are present) is indeed for 3/4 filling. Because 
of Eq.flTOD the 3 electrons per lattice site are distributed in the ground state in such a 
way, that on every site we have always at least one / electron present. As a consequence, 
H = H + U = [Hq + U J2i,a Ma — U L] + U P' has the ground state wave function given 
by Eqs.(P,pH|), provided the energy of the / level in H is renormalized as E't = Ef + U 
and Eg is shifted down with U L. We underline that in the U ^ limit, the described | ipg ) 
becomes only a small contribution from the linear combination of wave functions that build 
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up the ground state. The physical properties of the interacting ground state are present 
only at U > 0, and the U 7^ state cannot be perturbatively obtained from the U — 
case. This confirms the general belief that a non-Fermi liquid emergence has to be a clear 



non-perturbative effect [21 



For U > 0, instead of H we have H = -B + rjN-UL + UP' and Eqs.(@ -0) remain 
all valid, excepting Eq.@, which becomes Ef = Ef + U = \a\ 2 + \(3\ 2 — I7I 2 — \S\ 2 . Because 
Eqs.(§-0) remain the same, the a, j3, 7, 5, 77, m values remain unaltered (together with e\, e 2 
in Eq.(H) below), we obtain H' Q = Eg + P + , where E u g = L 77 ( 1 - 2 |m| 2 ) - U L. 
The term P + = P + U P' is a positive semidefinite operator, whose minimum eigenvalue 
is given by | ip g ) since P + \ip g ) = 0. The conditions that enable the construction of the 
solution for H change from Eq.(§) to 

Ef = Ef + U = 77 ( 1 — |m| 2 ) , (11) 

r) and Vj 2 remaining as given in Eq.(|8|) and U > 0. We denote the manifold defined by 
Eq.flrj by£> p @. 



In fact I ip g ) can be written in k space as well. Fourier transforming and taking e\ = 
a* + (3* exp(ifc), e 2 = a* m + (3* m* exp(ik), = £ fc e ikr > F fc (a) , D kj(T = e x c\ a + 
e 2 Fl 3) = E CT a CT with if } = 0, F^ = - F^ , and M = 

J2p(— l) p exp[i(ri h lx + r 2 k i2 + ... + r^/c^ )] where J2p denotes a sum over all permutations of 
(1, 2, ... L) to (ii, 12, ■■■ih) and p represents the number of pair-permutations in a given P, we 
find \ip g ) = M 3 Uk (ilLi F i a) ) 1°) • With A = M 3 we finally obtain \ip g ) = AjJk 4|0), 
where 

dk = e\ 4 ;T c{ ;l F® + e x e 2 T £ a CT , (12) 

and (ipglipg) = \A\ 2 H k |e 3 | 2 [ |ei| 4 + |ei| 2 |e 2 | 2 ] . The e 3 constant coefficient is chosen to 



preserve the normalization to unity |25[ and 1 63 1 2 = J2 a \ a cr\ 2 - 

To have an insight about the physical behavior of the system in the ground state, let us 
diagonalize first Ho in momentum space. Denoting by H k = ( V k c ka . f k>(T + h.c. ), in the 
thermodynamic limit we find 



where e% 



Ho = 
L 

2 t c cos k, e{ 



E 



2 ^ dfc 



o 2vr 
Ef + 2tf cos fc, Vk 



E 4 &L V* + H > 

b=c,f 



(13) 



2 V\ sin A; + z Vq. Introducing the row 



vector Wl = ( c ka , ) and the ( 2x2 ) matrix R with components R(i,i) = e c k , R(x,2) — 
Vk, R(2,i) — V k * , R(2,2) — e ki we ma y write the integrand in Eq.flTBD as W\ R Wk- The 
diagonalization in k space reduces to the secular equation written for R. Two bands arise, 
that via Eq.@ becomes 



E 



(i) 



77 2 



+ t > 



E 



(2) 



2 a/w 2 + tl + 2 ( t f + t c ) cos k 



(14) 



where we have as presented in Eq.(H), sign(tf) = —sign(t c ). The system is at 3/4 filling, so 
the lower band E^ is completely filled, and the upper band E k \ which is completely flat, 
is half filled. There is no hybridization between these two bands and taking into account 
their filling, the ground state energy E g is re-obtained. In the presence of U > 0, in the 
ground state, because P 1 \ ip g ) = 0, the effective Hamiltonian is in fact H' Q / L which differs 
from Hq/L in that it has a renormalized E't = Ef + U, and its energy scale is shifted down 
with U. Effectuating the calculations for the band structure as presented for Eq.(|l4]) but 
using instead of Ef the Ef value, we re-obtain (shifted down with U) for the ground state 
the structure presented in Eq. (|14|) . So we have a well defined Fermi energy, positioned at 



E 



(i) 



constant, the Fermi momentum being undefinable. The system has a natural built 



in gap (the minimum and nonzero distance between the E^ and E^) that will be visible 



in the physical properties at T ^ |26j 



To understand the physical behavior of the system all ground state expectation values 
at U > relevant for our study can be expressed from Eq.([TJ). Introducing Z^ 1 = 



l e i| 4 + l e i| 2 l e 2| 2 ) l e 3| 2 , we have 



\.o fk,a ) 



) = 


Kl 2 


leil 4 


+ | 


) = 


Kl 2 


|ei| 2 


M 


) = 




2 |ei 


2 * 
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Zk , 



(15) 



It is seen that = J2a ( ( c l l(T c k,a ) + ( ft yCr fk,a ) ) = 3 i.e. the total momentum 
distribution function is uniform in k space, so the system is a non-Fermi liquid. From 



Eq.(15) all individual contributions in can be expressed. We obtain functions of k that 
are continuous together with their derivatives of any order in the whole momentum space. 
For example, n c k ^ = ( c£+ Ck,i ) = 1 — A (rj + 2 t c cos k)(fj + 2 t c cos A;) -1 , where the 
constants are A = \a^\ 2 \m\ 2 / [|e 3 | 2 (|m| 2 — 1)] and fj = i] (|m| 2 + l)/(|m| 2 — 1) > 2t c . 
The movement of particles is allowed by \i> g ), that requires ,,at least one / electron on 
every site", i.e. allows nonzero hopping matrix elements. Indeed, starting from Eq.(|l5|) 
the ground state expectation values of all contributions to H can be expressed. With 
I / L = J 2 " dk[2n{\e 1 \ 2 + |e 2 | 2 ) J^ 1 > 0, we find 

(T c ) , ,9 (ff) , ,. (E f ) 



A 1 " 1 ' A 1 " 1 ' A 

(U) = U (A 2 - L), (V ) = -2V 2 I 



2r] 


m 


2 (1 


+ 


\m 


2 )L 


(1 " 


m\ 


2\2 



Vi) = ~ ?i r-U + ^3 , (16) 



where A x = [L-rj(l+ |m| 2 )/]/(l - |m| 2 ) 2 , A 2 = [( 1 - 2 |m| 2 ) L + 2rj \m\ 2 1] / ( 1 - |m| 2 ), 
and A 3 = {2t] 2 [(1 + \m\ 2 ) / (1 — \m\ 2 )} 2 — 8 V 2 } I . Summing up all contributions in Eq. (^6|) 
we re-obtain exactly the ground state energy Eg . We observe from Eq.(^) that E^ cannot 
be expressed as a sum of on-site contributions. With |m| 2 = \m\ 2 ( 1 + |m| 2 ) / ( 1 — \m\ 2 ) 2 , 
J = rj ( 1 + |m| 2 ) I — L, (H loc ) = ( E f ) + ( U ) + ( V ), we have 



{E f ) + (U) + (H h ) = Eg +2 V \m\ 2 J > E v g , 

(H loc ) = E" + 8 V, 2 I > Eg . (17) 

Here J = L [ 1 - r~ 2 ]"V2 — L > and r = ( rj / 2 t c ) [ ( \m\ 2 + 1 )/( |m| 2 - 1 ) ] 
. From Eq.([HJ) it can be seen that the system is not localized because (Hi oc ) > Eg 7 , i.e. 
the sum of the ground-state expectation values of all on-site (localized) terms from the 
Hamiltonian is greater than the ground-state energy itself. On the other hand, from Eq.(|T6D 
J2b=c,f (Tb) = — 2 rj \m\ 2 J < . Furthermore, as can be seen from Eq. (|TE|, |TTD adding 
the non-local part of the hybridization to the kinetic energy contributions (T c ) + (Tf) we 
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obtain also a negative number (H mov ) = (T c ) + (Tj) + (Vi) < 0, [p3| . So taken together, the 
sum of all ground-state expectation values connected to the movement of particles within 
the system is negative. As a consequence, similar to the 2D case with flat upper band and 
imaginary m parameter |24], in the situation described here, the system is maintaining its 
itinerant character instead of becoming insulating, because the (H mov ) (i-e. the movement 
of particles within the system) decreases the total energy As a consequence, the system is 
in a normal phase, but without having any non-regularity in n k at any point in k space (even 
in its derivatives with respect to k). From Eq. (|T0D , the coefficient remains undetermined, 



the ground state has a large spin degeneracy, so it is paramagnetic (see also Ref. f25|j). 

From the point of view of the described properties is important to analyze the system 
also around D p . Concerning these aspects, first of all it can be seen that the presented 
state emerges also in the presence of doping. Taking N = 3L + n r , where 1 < n r < L 
represents the additional electrons introduced into the system above 3/4 filling, we define 
^ (4) = E W A {k} nr /f with = e T (/+ T + e«*»c£ T ) + e^Uti + e^4j , where A {k} 
represents the coefficient of a given {k} combination of n r different and ordered ki taken 
from the L possible k- values. The ground state wave function becomes under doping \ip g d) = 
n^nLi^)^! ); the character of the ph ase described remaining unchanged since the 
first three operatorial components from \ip g d) remain as in Eq.(^). As can be seen, the state 
we are describing persists also above half-filling concentration of particles within the upper 
band. Concerning the N dependence, we further mention that the system is a Fermi-liquid 
for N < 2L. 

Taking into account that the flat-band feature is intimately connected only to the 
presence of the B term in the Hamiltonian (B collecting all k dependences), we can 
maintain the described flat-band features together with the non-Fermi liquid properties 
taking into consideration some local or global contributions to H that keep the system 
itinerant (otherwise the system becomes insulating). For example, we can leave the D p 
manifold taking into account arbitrary and independent deviations 5Ef and 5U from 



9 



the Ef and U values that satisfies Eg. (|TTD . Denoting by 5X = SEf + SU, considering 
U' = U + 5U > 0, the new Hamiltonian becomes H> = H + 5XJ2iaMa- Adding and 
subtracting 5XN C , a new positive semidefinite operator emerges (we consider here 5X > 0) 
Px = 5XEi, ff (l — c ?> c ?,o-), an d we have H' = (P + + Px) + E' g , where taking also doping 
into account, E' g /L = r][l — 2\m\ 2 + (n r /L)(l + 5X/i])] — U' . The new ground state wave 
function is obtained as P c \ip' gd ), where in \ip' gd ) the term introduces only /-electrons, and 
P c = riio-^icr maximize the number of c-electrons |27|- Under these conditions, at n r = 0, 



the obtained ground state wave function is a Mott insulator since it sets rigorously three 
electrons on every site of the lattice (similar result is obtained at 5X = n r = and real 
hybridization coupling constants). As a consequence, the described non- Fermi liquid state 
emerges in parameter space in the close vicinity of an insulating phase. In the presence of 
doping, at 5X > the on-site expectation value {Hi oc ) = (r]—\m\ 2 + 5X)(L+n r ) — U'L > E g , 
the system is itinerant, and since the flat-band feature is present, it is again non- Fermi liquid. 

To further emphasize the itinerant character of the ground-state described here we are 
presenting a localized ground-state obtained at 3/4 filling in another region of the phase 
diagram, but the same model, corresponding also to a completely flat upper band 



iv*> = ni«T + + + + ^4t + 5 4i)ii°) > ( 18 ) 

i 

where cti, A>7i)^ are constant numbers. This ground-state contains the same number of 
particles (i.e. 3) on every lattice site and gives |3S 

(t) = (f» = (Vi) = , (19) 

(H loc ) = Eg t , (20) 

where, for the localized case, the ground-state energy per site is E^J L = —(V Vi/tf)(l — 
2t'j/Vi) — U . A comparison of the two ground-states from Eq.(|9]) and Eq.(|I8D shows that 
sets rigorously three electrons on each site of the lattice, while \ip g ), besides sites with 
three electrons, contains sites with four and two electrons as well. As a consequence, any 
hopping like Hamiltonian term (as T c , Tj or the ,, band-change hopping" Vi) applied to 
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gives a state orthogonal to from where Eq.(^) arises. This is no more true in case of 
\i/j g ), from where the main difference between the two ground-states emerges. From the point 
of view of the Kubo formula, we mention, that starting from it we have (see for example 
p9| ) 1/m* ~ —(Hkin) and the real part of the conductivity in u — > limit o'{uo) ~ 1/m*, 
where m* is the effective mass, and Hkm is the global kinetic energy (in our case build up 
from T c , Tf and V% which give contribution in the current operator). As a consequence, we 
have a localized state only if m* = oo, i.e. in case where Eq.(|I"9D or Eq. (|20|) hold. 

In conclusion, in this Letter a new non-Fermi liquid normal phase has been presented, 
which is not intimately connected to one spatial dimensions. 
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